Compaction of anisotropic granular materials 
experiments and simulations 
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We present both experimental and numerical investigations of compaction in granular materials 
composed of rods. As a function of the aspect ratio of the particles, we have observed large vari- 
ations of the asymptotic packing volume fraction in vertical tubes. The relevant parameter is the 
ratio between the rod length £ and the tube diameter D. Even the compaction dynamics remains 
unchanged for various particle lengths, a 3d/2d phase transition for grain orientations is observed 
for £/D — 1. A toy model for the compaction of needles on a lattice is also proposed. This toy 
model gives a complementary view of our experimental results and leads to behaviors similar to 
experimental ones. 

PACS numbers: 05.70. Jk,45.70.Cc 



INTRODUCTION 

Granular matter has been the subject of numerous 
studies since the last decade Q, S 0- Indeed, most 
of the industrial products are processed, transported and 
stocked in a granular state. The packing density of those 
granular materials becomes therefore a relevant param- 
eter for a broad range of applications. The best way to 
reduce the costs for the manipulation of such granular 
materials is to increase the packing density p. This can 
be achieved by tapping or vibrating the vessel containing 
the grains. 

Various experimental studies 0,0 have underlined the 
fact that the dynamics of compaction is a complex prob- 
lem. The compaction dynamics is indeed characterized 
by a slow dynamics 0. Different laws have been pro- 
posed for the volume fraction evolution of a granular 
material as a function of the number n of taps. Among 
others, one has proposed the inverse logarithmic law 



p{n) = Pc 



Ap 



l + B ln(l + 2) 



(1) 



where the parameters p^ and Ap are respectively the 
asymptotic volume fraction and the maximum variation 
of the volume fraction. The dimensionless parameter B 
depends on the acceleration during each tap and t is 
the relaxation time of the reorganization process. This 
inverse logarithmic law of n was obtained in numerical 
models for compaction like the Tetris model and could 
also be derived from theoretical arguments Q. 

The great majority of earlier experiments have consid- 
ered the compaction dynamics of spherical grains H . 
Only a few papers discuss the problem of anisotropic par- 
ticles 0, ^3 ■ It has been reported that the compaction 
dynamics exhibits different regimes associated to respec- 
tively grain translations an d g rain rotations. 

It has been also reported [lillldl l that the random pack- 
ing density of spherocylinders or rods decreases when the 
aspect ratio £/d of these objects increases. 



In this paper, we present an experimental study of 
compaction dynamics for cylindrical particles in a ver- 
tical tubes (Figure 0). We investigate the volume frac- 
tion as a function of the aspect ratio of the grains £/d 
and as function of the ratio between the grain length and 
the tube diameter £/D. We will see that according to 
the ratio t/D, restrictions on grain orientations appear. 
These restrictions have large influence on the asymptotic 
volume fraction. 

A toy model for compaction is also proposed for the 
case of anisotropic grains. 



EXPERIMENTAL STUDY 
Experimental set-up 

Our experimental setup for the study of compaction 
is different from usual ones. Indeed, an oscillating sys- 
tem is generally used to produce the compaction process. 
In our experiment, an electromagnetic hammer, which is 
controlled by a computer via an interface, is placed below 
a vertical tube for tapping. Thus, we act on the system 
with a series of intense and short taps. The intensity, the 
number and the frequency of the taps can be controlled. 
A CCD camera records the tube during the whole exper- 
iment. Therefore, we can measure the evolution of the 
height h of the granular/air interface by image analysis. 
A ruler is placed along the tube in order to calibrate the 
measurements. The picture analysis for the volume frac- 
tion measurements is made by an algorithm especially 
developed for this task. On the pictures, grains appear 
in dark and the rest of the system appears in bright. An 
average of the pixel values is performed on the vertical 
of the picture. We obtain an abrupt variation of the av- 
erage at the position of the granular/air interface. The 
maximum of the first derivative gives the position of the 
interface. The volume fraction p is estimated through 
the value of h, the density of grains and the weight of the 
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entire column. One should note that in our experiment, 
we are mainly interested on the volume fraction averaged 
over the entire granular column. A sketch of our setup is 
given in Figure ^ 




FIG. 1: Sketch of our experimental setup. The granular ma- 
terial is placed in a vertical glass tube. Taps are generated 
by an electromagnetic hammer controlled by a computer. A 
camera records the granular/air interface after each tap. 



The granular material we used is made of identical 
cylinders. The diameter d of these cylinders is 1.2 mm 
and their length varies between 2 mm and 26 mm. Their 
aspect ratio i/d varies between 1.6 and 21. The cylinders 
are placed in a vertical glass tube. The internal diameter 
D of the glass tube can vary between 8 mm and 24 mm. 

For the preparation of the loose packing, grains are 
poured in the glass tube with a funnel before each run. 
We have also tested a decompaction by inflating air from 
the bottom of the tube. This does not change our results 
obtained with the rain method. 

In order to obtain a saturation of the volume fraction 
when tapping, a minimum of 2000 taps was applied be- 
fore stopping each experiment. 

Before investigating the compaction dynamics, let us 
present some additional information about tapping. We 
have measured the acceleration experienced by the bot- 
tom of the glass tube after each tap. The analysis of a 
single tap with an ultrafast camera gives a typical ac- 
celeration of 12 g during a short period of 2 ms. The 
amplitude of the vertical motion is approximately 0.25 
mm. We have also measured the relaxation time of the 
system after a tap (the time after which there is no move- 
ment in the tube). We obtained at most 0.12 s. In all our 
experiments, the successive taps are always separated by 
at least 0.2 seconds. This time interval is long enough to 
avoid any overlap between relaxation mechanisms pro- 
voked by two successive taps. 



Compaction dynamics 

We have analysed the compaction for different tube di- 
ameters D and for different lengths I of the rods. The 
figure presents typical curves of compaction for three 
different lengths. Each curve is well fitted by the inverse 
logarithmic law of EqQJ. Curves obtained in other glass 
tubes present similar behaviors. We do not observe dif- 
ferent stages in the compaction process as in [lll |. In fact, 
several stages have been observed for a vertical accelera- 
tion less than T = 7.5. In our experiment, the dimension- 
less acceleration is much larger (r = 12) as underlined 
here-before. 

0.66 I ■ — ■ — ■ — n 
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FIG. 2: Experimental compaction curves : volume fraction p 
as a function of the tap number n. The curves are fits using 
Eq. Q. Three different cases are illustrated : I = 4 mm 
(squares), t = 10 mm (circles) and t = 20 mm (triangles) for 
D = 10 mm. 

All experimental compaction curves have been fitted 
with the inverse logarithmic law, i.e. with Eq.(l). From 
the fits, a relevant parameter is the asymptotic volume 
fraction p^. The figure [3] shows the asymptotic volume 
fraction p^ as a function of the ratio i/D for different 
tube diameters. The asymptotic volume fraction presents 
large variations when the ratio i/D is changed. In all ex- 
periments, we observe a minimum of p^ when the cylin- 
der length is equal to the tube diameter, i.e. when I = D. 
From these results, one could ask whether the compaction 
dynamics changes or not when varying i/D. 

While the asymptotic volume fraction of the pack- 
ing does not vary much with D for both extreme cases 
i/D >> 1 and i/D « 1, the situation is quite dif- 
ferent for the particular case I = D. In that case, the 
asymptotic volume fraction p^ decreases when the tube 
diameter D increases. Figure 01 is a plot of p™ m as a 
function of I for the particular case i — D. The line in 
Figure 01 is only a guide for the eyes. 
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FIG. 3: Asymptotic volume fraction p^ as a function of the 
ratio l/D for three different tube diameters : D= 8 mm 
(squares), D = 10 mm (circles) and D = 12 mm (triangles). 
Error bars are indicated. 
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FIG. 4: Asymptotic volume fraction p™" as a function of I 
for the particular case i — D. Error bars are indicated. The 
curve is only a guide for the eyes 



Discussion 

As for common liquid cristals, the rods have two pa- 
rameters for describing a possible long-range order : po- 
sition and orientation. Each of them plays a role in 
the packing fraction. Since random positions are usu- 
ally observed even for spherical grains, the largest source 
of heterogeneities in our experiment comes from grain 
orientations. The orientation disorder is reduced when 
£ > D. One should also note that the vertical configu- 
rations which are favored by the geometrical constraints 
do not correspond to energy minimization. Indeed, the 
potential energy for a grain is minimized when this grain 



FIG. 5: Pictures of the tube after 2000 taps and for a tube 
diameter D — 10 mm. Three different cases are illustrated: 
(a) I = 4 mm, (b) I = 10 mm, (c) I = 20 mm. 

Figure presents typical pictures of the packing for 
various ratios l/D. The pictures were taken after sev- 
eral taps. From the observation of the left picture shown 
in Figure one understands that all grain orientations 
are possible for I < D. In that case, the ideal packing 
could be realized in various directions. For £ > D, hori- 
zontal orientations are forbidden. The grains are ideally 
arranged along the tube direction, i.e. along the vertical 
direction. This global ordering of the grains is well seen 
on the right picture of Figure 03 When £ = D, a few 
grains are placed in the horizontal direction. Large voids 
can be seen near horizontal grains. This case is shown in 
the central picture of Figure |5] 

The large variation of the packing fraction, which is ob- 
served when the needle length £ is equivalent to the tube 
diameter D, should be attributed to a 3d-2d transition. 
Let us discuss some particular cases. For fCfl, the situ- 
ation tends to a packing of isotropic objects. The 3d ideal 
volume fraction of identical spheres is /? = 7r/3\/2«0.74 
while the random close packed limit is p ps 0.58, i.e. a 
value close to what we have measured for the smaller 
grains. In the case of long needles (£ S> D), all needles 
have a vertical orientation. The ideal volume fraction is 
given by the fraction of disks arranged along a hexag- 
onal packing in the plane crossing the glass tube, i.e. 
p = ttV3/6 ps 0.91. Of course, this ideal ordered pack- 
ing is never realized. The random case corresponds to a 
packing fraction p ps 0.82 which is the upper limit of our 
experimental measurements. 

From the above discussion on both extreme cases, one 
expects a global increase of the volume fraction from 0.58 
(random sphere packing) to 0.82 (random disk packing) 
when the ratio £/D increases. How to explain that vol- 
ume fraction presents a minimum for £ ps £)? This mini- 
mum is the consequence of the competition between two 
phenomenons. First, the packing density decreases when 
the aspect ratio of grain l/d increases [I3j,llj. Secondly, 
we observe an arrangement of grains along the tube di- 
rection when l/D > 1. This arrangement leads to an 
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increase of the packing density. The first phenomenon 
explains the decrease of poo for £/D situated between 0.2 
and 1. The second phenomenon explains the increase 
of Poo for £/D situated between 1 and 2. Furthermore, 
when £ sa D, some horizontal grains can block the verti- 
cal motion of the other grains and form jams. 

In Figure 0] the minimum value of p^ is found to de- 
crease with £ (=D). This decrease is explained by the 
increase of the aspect ratio £/d and by the presence of 
jams when a grain is horizontal. 



SIMULATIONS 



gravity when randomly selecting new orientations for the 
needles. Indeed, the effect of gravity becomes negligible 
when p becomes large (after a few taps). 



Toy model 




In order to reproduce our experimental results, we pro- 
pose a toy model on a 2d lattice. Grains are needles hav- 
ing different orientations. They are placed on a square 
lattice. Each grain has 8 different possible orientations 
which arc illustrated in Figure 7. Lattice boundaries are 
closed. The length £ of the needles is a parameter to be 
investigated. The width W of the lattice is the second 
parameter of the toy model. 



FIG. 6: Eight possible orientations of a needle with £ — 9 on 
a square lattice. 

The first step of the model consists in filling the lat- 
tice with anisotropic grains. The first configuration is 
obtained by the rain method. Grains are dropped se- 
quentially on the top of the lattice. Both grain position 
and grain orientation are chosen randomly. The grains 
follow a vertical motion until they touch another particle 
or the lattice edge. 

Once the lattice is filled, the model consists in apply- 
ing a scries of taps. The procedure to simulate a single 
tap is the following. We randomly choose some grains 
as many times as the number of grains in the system. 
Thus, every grain in the system has a chance to move. 
For a selected grain, we have the choice between two mo- 
tions : rotation or translation. The selected motion for 
that grain has also two possible directions. Each motion 
and each direction is chosen randomly with a probability 
0.5. The operation is executed only if it does not lead 
to some overlap of the grains. The volume fraction p of 
the packing is measured after every tap. It should be 
also underlined that this toy model does not consider the 



FIG. 7: Typical simulated packings in a 2d system. Three 
different needles {I = 13,41,77) are illustrated for the same 
lattice (W = 41). Different heterogeneities from disorder to 
long-range order can be seen. Large voids are obtained when 
£ = W. 



Numerical results 

In the simulations, we have observed the evolution of 
the volume fraction p function of the tap number 
n. Figure |S] presents typical compaction curves obtained 
numerically for 3 different situations. An inverse loga- 
rithmic dynamics has been found in all cases. From the 
fits using Eq.(l), we have obtained the evolution of the 
saturation density for various ratios £/W and for dif- 
ferent lattices. The results are displayed in Figure El As 
for experiment the volume fraction presents a minimum 
when £ = W. 

Moreover, we have numerically investigated the partic- 
ular case £ = W. Figure ITUI presents the volume fraction 
Poo as a function of W in a semi- log plot. The volume 
fraction is seen to decrease exponentially with the lattice 
width W . One has 



W 

Pr - A exp( — — ) 



(2) 



with pn = 0.125 ± 0.009 which is a residual density for 
a jammed phase and L = 21.1±2.3 is a characteristic 
length of particle. 



Discussion 

Our toy model reproduces qualitatively our experimen- 
tal results. Indeed, the compaction curves look similar 
to the experimental curve (see figure |SJ) . A minimum 
of the volume fraction is also observed when the needle 
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FIG. 8: Simulated compaction curves: volume fraction p as 
a function of the tap number n. The curves are fits using 
the inverse logarithmic law Eq. Q. Three different cases 
are illustrated : £ = 9 (squares), £ = 21 (circles) and I = 41 
(triangles). All data are obtained for W = 21. 



0.9 p 

0.8 - 

0.7 - 

0.6 - 

0.5 - 

0.4 - 

0.3 - 

0.2 - 

0.1 - 



-h: 



0.5 



1 

l/D 



1.5 



FIG. 10: Asymptotic volume fraction p^ as a function of £ 
for the particular case £ — W. Error bars are indicated. The 
curve is only a guide for the eyes 



large voids. In the third case, a large majority of the 
grains have a vertical orientation and packs ideally along 
the lattice. A few diagonal grains creates voids in the 
lattice. 



CONCLUSION 

In this paper, we present an experimental study 
of compaction in granular materials composed of 
anisotropic particles. As a function of the aspect ra- 
tio of the particles, we have observed large variations of 
the asymptotic packing volume fraction in vertical tubes. 
We have observed a transition 3d/2d when the needle 
length becomes equal to the width of the container. A 
toy model is proposed and reproduces numerically what 
we observed. 



FIG. 9: Asymptotic volume fraction poo as a function of the 
ratio £/D for different lattices: W=13 pixels (squares) and 
W=41 pixels (circles). Error bars are indicated. 



length I becomes equivalent to the width W of the lat- 
tice. The exponential scaling of the minimum (figure rTU|) 
is unexpected and needs further theoretical studies. 

The large variation of the packing fraction, which is 
observed when the needle length £ is equivalent to the 
width W of the lattice, should be attributed to a 2d- 
ld transition. Figure presents typical packings on the 
same lattice. Three different lengths I are illustrated: 
£ < W, I = W and I > W . In the former case, the grains 
have various orientations. Small voids can be seen. In 
the second case, a few grains have horizontal orientations 
blocking the lattice. This jammed situation could create 
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